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PERSISTENCE OF SADDLE BEHAVIOR IN THE NONSMOOTH LIMIT OF 

SMOOTH DYNAMICAL SYSTEMS 

JULIE LEIFELD 1 , KAITLIN HILL 2 , AND ANDREW ROBERTS 3 


Abstract. Models such as those involving abrupt changes in the Earth’s reflectivity due to ice melt 
and formation often use nonlinear terms (e.g., hyperbolic tangent) to model the transition between 
two states. For various reasons, these models are often approximated by “simplified” discontinuous 
piecewise-linear systems that are obtained by letting the length of the transition region limit to 
zero. The smooth versions of these models may exhibit equilibrium solutions that are destroyed 
as parameter changes transition the models from smooth to nonsmooth. Using one such model 
as motivation, we explore the persistence of local behavior around saddle equilibria under these 
transitions. We find sufficient conditions under which smooth models with saddle equilibria can 
become nonsmooth models with “zombie saddle” behavior, and analogues of stable and unstable 
manifolds persist. 


I 


1. Introduction 

One of the fundamental concepts in the analysis of smooth dynamical systems is the idea of 
local behavior. However, varying a parameter can transform a smooth system into one that is 
nonsmooth, and nonsmooth systems often have local behavior that is not observed in smooth 
systems. Indeed, in many nonsmooth systems, even the standard uniqueness theorem does not 
apply. In this paper we analyze persistence of saddle-like behavior in nonsmooth systems. We 
conclude that if a saddle equilibrium disappears as a result of the system becoming nonsnrooth, 
and not in the case of a standard smooth bifurcation, saddle behavior will in some sense remain 
in the nonsnrooth system. We show that the saddle equilibrium in the smooth system limits to a 
nonsmooth pseudoequilibrium, which is an equilibrium of a Filippov flow. This pseudoequilibrium 
is accompanied by stable and unstable manifolds that are analogous to their smooth counterparts. 
We call the new pseudoequilibrium point a zombie saddle , and say the system has zombie saddle 
behavior. 
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The system given by 


/ p(x) M x ) > 0 
1 q(x) M x ) < 0 , 


where p(x) and q(x) are smooth functions, is an illustrative example of a nonsmooth system. Notice 
that the dynamics on the manifold where /i(x) = 0 remain undefined. This manifold where h(x) = 0 
is called the discontinuity boundary , or splitting manifold , and the system is said to be a piecewise- 
smooth system. We are particularly interested in relating the dynamics of a smooth system to 
those of a nonsmooth approximation to the system obtained through some limiting process. 

Conceptual climate models have been formulated using nonsmooth systems dating back at least 
to Stommel’s thermohaline circulation model in 1961. Over the last decade, the use of nonsmoth 
systems has become increasingly popular, with [B ia ia lama Hues! all discussing conceptual climate 
models that are either nonsmooth or exhibit nonsmooth behavior in some limit. In particular, a 
nonsmooth approximation is often used for the basic energy balance equation 


(1) — = Q(1 -a(T))-(A + BT). 

The albedo function a(T) measures the reflectivity of the Earth’s surface. It is generally assumed 
to transition from a high albedo at low temperatures due to large quantities of ice to a low albedo 
at high temperatures. The transition is often modeled using a hyperbolic tangent: 


( 2 ) 


a{T) = 


oti + or 


Oid QLq, 


tanh 


T-T 0 

D 


D measures how rapidly the transition from high to low albedo occurs. In the limit D —> 0, the 
transition is instantaneous, and a is discontinuous. In (fTj), the dependence on greenhouse gas 
concentration is built into the parameter A, with smaller A indicating more greenhouse gases. 
However, especially on the time scales of paleoclimatic phenomena like Snowball Earth, A should 
not be treated as a constant Gang. Introducing A as a (possibly slow) dynamic variable gives 
the system 


f = Q(l-a(T))-(A + BT) 
A = —mA + nT + C, 


where m,n > 0 [!J. For certain parameter values the system has three equilibria (see Figure GK a)), 
however as D —> 0 the middle equilibrium disappears, as illustrated in Figure (D(b) . In the case 
when the smooth system is bistable and the middle equilibrium is a saddle, the stable manifold to 
the saddle forms the boundary of the basin of attraction of the two stable equilibria, corresponding 
to stable climates. 

In the context of nonsmooth systems, there has been significant interest in the behavior of 
systems as they transition from smooth to nonsmooth. For a more in-depth treatment of the recent 
developments in the dynamics of nonsmooth systems, we refer the reader to the surveys by Colombo, 
et al. [Jj and Makarenkov, et al. [18]. As a system with multiple time scales transitions from smooth 
to nonsmooth, a technique called pinching has been proposed as a method of preserving dynamics 
of the smooth system in the nonsmooth limit [20l l8j. On the other hand, techniques such as blow 
up and regularization use geometric singular perturbation theory to approximate the behavior in 
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Figure 1. Nullclines of ([3]) when m = n = 1, Q = 300, cu = 0.8, a w = 0.2, 
To = 290, B = 4, C = —1300. (a) shows the case for D = 5, and (b) shows the limit 
as D —> 0 with the splitting manifold T = Tq (dashed, red line). Notice the middle 
equilibrium from (a) is no longer an equilibrium in (b). 


the discontinuity boundary of a nonsmooth system as the system is smoothed [211116] . This process 
is similar to the blow up method used by Jeffrey m, except in |15| more complicated dynamics 
are allowed to occur in the blown-up splitting manifold. 

Nonsmooth systems in which an equilibrium collides with a discontinuity boundary as a param¬ 
eter varies are said to undergo a border-collision bifurcation mm [Ej. However, the interaction 
between an equilibrium colliding with a discontinuity boundary as the system itself is limiting to 
a nonsmooth limit is not well-understood. Smooth systems that limit to a nonsmooth variant 
as some parameter varies have held recent interest, notably in the context of climate dynamics 
[23 , 3| [10, [23]. In this paper we specifically study systems with saddle equilibrium points that 
disappear in a traditional sense as the system becomes nonsmooth. By this we mean that the 
equilibrium limits to the nonsmooth boundary as a —> 0, and hence is not a place where the vector 
field defined in the system is 0. Before completing the analysis, we first define some syntax in 
Section [2j In Section [3] we state and prove the main results of the paper. We consider three cases of 
systems along with their nonsmooth limits: a system that has no sliding region in the nonsmooth 
limit, a system with a repelling sliding region, and a system with an attracting sliding region. In 
all cases we prove the existence of saddle equilibria in the systems throughout the limiting process. 
Additionally, we show show that the saddle behavior persists in the limit. We call these persistent 
pseudoequilibria zombie saddles, because they are the remnants of “real” saddle equilibria from 
the smooth systems. Then in Section 0] we provide examples of systems which exhibit this type of 
behavior. We conclude in Section [5] with a discussion of our results. 
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2. Nonsmooth Analysis Techniques 


Our analysis uses the standard syntax of nonsmooth systems. For simplicity we will assume 
there is a unique C 1 splitting manifold of codimension 1. Then, we can define a scalar function 
/i(x) : — > R, with h(x) = 0 along the manifold. In the regions where the system is smooth, we 

let / + (x) be the vector field where /i(x) > 0, and define / _ (x) similarly. Let A be a step function, 
defined 

. _ j 1 h(x) > 0 
| 0 /i(x) < 0. 

Then, as in £2], the entire system can be written in the form of a convex combination of the vector 
fields 


( 4 ) 


x = f(x, A) 


/ /+(x) h(x) > 0 
1 / (x) h(x) < 0 

A/ + (x) + (1 - A)/“(x). 


A sliding solution exists along the nonsmooth boundary if the following system can be solved for 
A £ [0,1]. 


( 5 ) 


S = f (x, A) • V/i(x) = 0 
/i(x) = 0 


In a Filippov system, solutions to Q are only possible in regions where the dot product / + -/~ < 0. 
Intuitively, one would expect sliding solutions here, as solutions should cross the splitting manifold 
when the vector fields point in the same direction. Therefore, sliding regions are defined as regions 
along the splitting manifold where f + -f~ < 0. The boundaries of the sliding region are then places 
where either vector field is tangent to the splitting manifold, or more formally, when ■ V/t = 0. 
A tangency in f + is called visible if a solution through the tangent point exists entirely in / + , and 
is invisible if the solution exists entirely in f~. Corresponding definitions exist for tangencies in 
f~. The stability of the sliding solution is determined by the sign of 


d 

d\ b 


along h(x) = 0, with 
indicating a stable sliding region. Alternatively. 


Tx s <° 




indicates that the sliding region is unstable. See Figure [2] for an illustration of attracting and 
repelling sliding regions. Recent work by Jeffrey has expanded on this method, by using the idea 
that convex combinations of / + and f~ are not necessary to the analysis [15]. For example, adding 
nonlinear terms of the parameter A, such as A(1 — A), which are zero away from the nonsmooth 
boundary, also gives a method of finding sliding solutions which are beyond the scope of the standard 
Filippov analysis. These solutions can exist in regions for which the vector field points in the same 
direction, and hence one would not expect sliding solutions. To avoid these difficulties, in this paper 
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we assume the nonsmooth functions can be written linearly with respect to A. In other words, we 
let g a (x) be a smooth, monotonically increasing function, and as a approaches 0, g a {x) approaches 
g 0 (x), with 


( 6 ) 


9o(x) 


g+(x ) x > 0 
g~(x) x < 0. 


It is possible to define go(0) based on a pointwise convergence of <y a (0). However, this precludes the 
possibility of any dynamic behavior along the splitting manifold which may be indicated by real 
behavior in the application system. Therefore, we instead use the nonsmooth analysis conventions 
previously described, and also assume 


9o (x) = Xg+(x) + (1 - A )g-(x). 

It is important to note the breakdown of uniqueness of solutions in sliding regions of nonsmooth 
systems. Several initial conditions can generate solutions that approach the nonsmooth boundary in 
either forward or backward time and intersect. In Figure [2j a), solutions are nonunique in backward 
time. On the other hand, in Figure [2jb) solutions are nonunique in forward time; that is, only one 
solution trajectory can enter the sliding region from either the right-hand or left-hand side, but a 
family of solutions leaves the sliding region from both sides. Hence, uniqueness in both forward and 
backward time is not to be expected in these systems. However, as we will see in Section [3j it may 
still be possible to choose unique trajectories which correspond to stable and unstable manifolds. 


3. Results 

In this section we prove the persistence of saddle equilibria in three distinct cases, when the 
sliding region is stable, unstable, or nonexistent. In each case, we prove the existence of saddle 
equilibria in the smooth systems, as well as saddle psuedoequilibria in the corresponding nonsmooth 
systems. We also prove that stable and unstable manifolds persist in some form. 

3.1. No sliding solutions. 


Theorem 3.1. Let 

9a{x) = g 

a/0, be a monotonically increasing, smooth function such that g a {x) —> go{x) as a —>• 0, with 

g+{x) x > 0 
9-{x) x<0, 

where g + and g- are smooth functions up to the boundary x = 0. Furthermore, assume go = 
Xg + (x) + (1 — A )g-{x), where 

1 x > 0 
0 x < 0. 


9o(x) = 


A = 


Finally, let g~(0) < /(0) < g+(0), where f is a smooth function. Then, when 0 < a < 1, the 
system 


( 7 ) 


x = y- f(x) 
y = 9a(x)-y, 


6 


J. LEIFELD, K. HILL, AND A. ROBERTS 


contains a saddle equilibrium point p a = (x a ,y a ), such that x a — > 0 as a —> 0. When a = 0, the 
nonsmooth system contains a pseudoequilibrium on the splitting manifold x = 0, and stable and 
unstable manifolds exist locally. 

Proof. First we look at the case o / 0. Equilibrium solutions are found by solving f(x) — g a (x) = 0. 
Consider an arbitrarily small interval around x = 0, (—6,5). Let d = /(0) — g+(0) < 0. Choose 
x\ £ (0, 5) such that 

|/(xi)-/(0)| < L 

ls+(o) -0+O»i)l < y 

It is possible to choose this x\ because of the continuity of / and g+. Then, choose a = ai / 0 
such that 

\g+(xi) ~ g a (xi)\ < y. 


1 


x 0 


-1 

-0.5 0 0.5 

x 

(A) 

Figure 2. Typical phase diagrams for nonsmooth systems with sliding regions, 
with (a) attracting sliding and (b) repelling sliding. In (a), trajectories may enter 
the sliding region at any point along the line between the dot and the square, while 
in (b) trajectories may only enter the sliding region at the dot from the left and at 
the square from the right. 


X 0 
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This small a exists because of the pointwise convergence of g a . We show that f(x\) — g a (x i) < 0. 
f(x i) - g a (x i) = /(0) - g+(0) + (g+{0) - g+(x i)) + (f(x i) - /(0)) + (g + (x i) - g a (x i)) 
<d+\d\ 

= 0 


A similar argument works for xi £ (—5,0). Therefore we conclude by the Intermediate Value 
Theorem that f(x) — g a (x) = 0 for some x £ (X 2 ,X\ ) C (—5,5). 

To see that this equilibrium is a saddle, it is convenient to make a change of variables, which is 
valid for all a/0. Let £ = x/a. Then, the system becomes 

i = -{y-m)) 

a 

y = g(0~v- 


The Jacobian of this system is 


with 


J = 


-/'«) 

g '(0 



det J = f'(aO - -g'(0- 

a 

Because g is monotone increasing, we know g'(0 > 0. Thus, 
sufficiently small o, det J < 0. 


it is immediately clear that for 


When a = 0, the system is discontinuous and the splitting manifold is the line x = 0. The 
boundaries of the Filippov sliding region occur where y = /(0). Given that / is a function, there is 
only one point at which this occurs, meaning that there is no linear sliding region. This point is a 
pseudoequilibrium of the Filippov flow, as the vector fields on either side of the splitting manifold 
are colinear and opposite here. It is easy to check that the the tangency is visibile on each side, as 
seen in Figure (3)( a). 

Furthermore, we can show the existence of trajectories through this tangency. First, we consider 
the system defined where x < 0, and the case /'(0) > 0. Let y > /(0). Then, because /(0) > g~(0), 
there is a region x £ (—5,0), y £ (/(0),/(0) + e) for which y < 0. Moreover, because f'(x) > 0, 
x > 0 in this region. This implies that there are no strictly invariant sets in the region, and therefore 
every point can be taken continuously to the left and upper edges of the boundary using the map 
defined by the flow. In particular, there is a trajectory in this set which maps to (0,/(0)), the 
tangency point, under the flow. This trajectory is an analogue of a stable manifold, along which 
the equilibrium can be reached in finite time. Because (0, /(0)) is a tangency of the vector field 
when x < 0, the trajectory remains in the x < 0 region, and the portion which approaches (0, /(0)) 
in backward time likewise acts as an analogue of an unstable manifold. If f'( 0) < 0, we can make 
a similar argument after reversing time. Then, there is a region x £ (—5,0), y £ (/(0),/(0) — e) 
for which x > 0 and y > 0. This implies the existence of in initial condition in this region which 
approaches the tangency (0,/(0)) in backward time, with respect to the original system. This 
trajectory will remain in the region x < 0 for the same reason, and acts as an analogue to stable 
and unstable manifolds. The same arguments can be made for the existence of a trajectory through 
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the tangency in the region x > 0. Thus in the full system, trajectories through the double tangency 
are not unique. Much like a saddle, there are two trajectories approaching the point, and two 
trajectories leaving. These trajectories function as stable and unstable manifolds to a saddle point, 
especially in their role as separatrices. Much like stable and unstable manifolds, their intersection 
consists of one point, however, in this case that point is reachable in finite time. See Figure [3jb) 
for an illustration of this idea, where the dark green boxes represent the boundaries of the flow in 
the regions discussed. □ 


2 


1 


is o 


-1 


-2 


Figure 3. (a) Example phase diagram of the nonsmooth limit of (|7|). with f(x) = 
2x + 0.4 and g a (x ) = tanh(~), so that go(x) = ±1. Blue curves represent the 
trajectories that have behavior similar to a stable manifold, and red curves exhibit 
behavior similar to an unstable manifold. The black points show the equilibria of 
the system. Dashed lines represent the nullclines of the system, (b) An illustration 
of the visible tangency using / and g as defined in (a). The dashed lines show the 
nullclines of the system, and the solid green lines give an illustration of the boxes 
considered in finding the tangency. 



3.2. A repelling sliding region. 

Theorem 3.1. Under the assumptions in Theorem 1, when 0 < a <C 1. the system 
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contains a saddle equilibrium point x a , such that x a —> 0 as a —> 0. When a = 0, the nonsmooth 
system contains a stable pseudoequilibrium in a repelling sliding flow on the splitting manifold x = 0, 
and stable and unstable manifolds exist locally. 


Proof. Again, equilibria are determined by the equation f(x)—g a (x) = 0. An identical calculation to 
that in the proof of Theorem 1 shows that an equilibrium solution exists in (—6, 6). The calculation 
will not be repeated here. Doing the same convenient coordinate change as before, we get a new 
system 


£ = 5 (g(£) ~ V) 
y = /«) - y- 


In this system the Jacobian is 


J = 


' Mo 

af(a£) 



with 

det J = /'«) - -g'(0- 

a 

Again we have that det J < 0, giving us a saddle equilibrium. 


The nonsmooth version of this system also shows saddle like behavior as g a (x) approaches a step 
function in the a = 0 limit. See Figure |4] for an example of a phase diagram of (J5|) . In this case, 
tangencies to the discontinuity boundary occur when y = g + ( 0) and y = g~(0), implying that the 
sliding region has nonzero length. Moreover, this sliding region is repelling, which means that any 
solution starting on the sliding region will leave arbitrarily. To see this, we will use some formalism. 
Recall that go(x) = A g+(x) + (1 — X)g-(x). We define h(x,y) to be a scalar function such that 
h > 0 when x > 0, and h < 0 when x < 0. Then, sliding solutions exist where f(x, A) • Vh = 0 
can be solved for some A £ [0,1] along x = 0. In this case, it is clear that sliding solutions exist 
on [<7_(0), <7_|_(0)], as asserted previously. Moreover, as in [15], stability of the sliding region can be 
determined by 

5(x,A) = M-Vh, 

where S < 0 implies an attracting sliding region, and S > 0 implies repelling. In this example, 
h = x, so 

f • V/i = \{g+{0) - y) + (1 - \){g-{0) - y). 

-^f ■X7h = g + (0)-g.(0) >0, 

indicating repelling sliding. 


We can also use this formulation to find the flow along the sliding region. Using (jSJ) to solve for 
A and plugging this solution into the original equations, we get an equation for the flow, 

V = /(0) - y. 

It is clear that this has one stable equilibrium solution, y = /(0). 
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Because the repelling sliding region has one stable equilibrium, the system has behavior analogous 
to a saddle. This saddle does not have unique trajectories acting as separatrices, as Section 13.11 
Instead, a family of trajectories will move along the sliding region, and each will leave the manifold 
at some point before it can reach the equilibrium. However, the sliding region itself might be 
considered analogous to a stable manifold. For the same reasons as described in section l3Jl there 
is a unique trajectory on each side of the discontinuity which emanates from the equilibrium on 
the sliding region. These two trajectories are analogous to an unstable manifold. □ 



x 


Figure 4. Example phase diagram of the nonsnrooth limit of (j8j), with f(x) = 
2x + 0.4 and g a (x ) = tanh(|), so that go(x) = ±1. Blue curves represent the tra¬ 
jectories that have behavior similar to a stable manifold, and red curves correspond 
to behavior similar to an unstable manifold. The points show the equilibria of the 
system, and the squares show where the sliding region begins/ends. Dashed black 
lines represent the nullclines of the system. 


3.3. An attracting sliding region. 

Theorem 3.1. Under the assumptions in Theorem 1, when 0 < a <C 1, the system 

x = y-g a {x ) 

1 j y = y-f(x), 

contains a saddle equilibrium point x a , with x a —> 0 as a —> 0. When a = 0, the nonsmooth system 
contains an unstable pseudoequilibrium in a stable sliding region on the splitting manifold x = 0, 
and stable and unstable manifolds exist locally. 
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Proof. The system in Theorem 3 is obtained by reversing time in the previous system, so we need 
not repeat the proof that a saddle equilibrium exists and limits to x=0. 

The nonsmooth system again has two tangencies along the discontinuity line. These occur 
at y = g+( 0) and y = g_(0). The sliding region in this case is attracting, because S(x, A) = 
g_(0) — g+(0) < 0. In this system, the sliding flow is given by the equation 

y = y- /(0) 

from which it is again immediately clear that the sliding flow has one unstable equilibrium. The 
equilibrium point is unstable, but has no unique unstable manifold. Instead, a family of trajectories 
will be attracted to the stable sliding region, and then slide away from the unstable equilibrium. 
The sliding region itself might again be considered as the analogue of the unstable manifold. Again, 
there will be unique solutions to the equations on either side of the discontinuity which hit the sliding 
region exactly at the equilibrium of the sliding flow. See Figure [5] These trajectories arrive at the 
equilibrium in finite time, and act as analogues to the stable manifold. So, again, the nonsmooth 
system displays saddle behavior, albeit slightly different from the saddles of the first two nonsmooth 
systems. □ 





X 


Figure 5. Example phase diagram of the nonsnrooth limit of ([9]), with f(x) = 
2x + 0.4 and g a (x ) = tanh(|), so that go{x) = ±1. Red curves represent the 
trajectories that have behavior similar to an unstable manifold, and blue curves 
correspond to a stable manifold. The points show the equilibria of the system, and 
the squares show where the sliding region begins/ends. Dashed black lines represent 
the nullclines of the system. 
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4. Examples 


4.1. A nonsmooth system that exhibits all three saddles. Finally, we find a nonsmooth 
system which displays all three types of saddle behavior. In this example we let g a (x) = tanh 
which limits to 

9o(x) = | 

The smooth system is as follows, where a > 0. 


x > 0 

x < 0. 


x = kg a (x)-x + y 
y = ag a (x) — y. 


-I + kg'Jxa) 1 
ag' a (Xa) “I ’ 


and det J = 1 — {k + a)g' a (x a ). This means that for a > 0 fixed and |fc| < a, the conditions given 
previously for g a {x) again guarantee that the equilibrium x a is a saddle point. 

The nonsmooth limit of the previous system is 


/ 10 n x = kgo(x) — x + y 

1 V = ag 0 (x) - y, 

where again 

go ( x ) = ( 1 x > = 2A - 1, 

\ -1 x < 0. ’ 

if A is defined as before. Tangencies to the discontinuity occur at y = —k( 2A — 1). See Figure 
[6] for illustrations of the phase diagram for this system for varied k. For this system, f • V/i = 
k( 2A — 1) — x — y. So, ^S(x, A) = 2k, and when k > 0, the sliding region is repelling. When k ^ 0, 
the sliding flow in this system is given by 

j=-(f + l) y . 

The equilibrium occurs at y = 0, and is stable if a > —k, which is true for any positive k. Similarly, 
when k < 0, the sliding region is attracting. In this case, the equilibrium is stable if a > |A;|, so in 
some region where \k\ is small, this is satisfied. 


Because of the division by k in the sliding flow, the case where k = 0 must be treated separately. 
In this case the system is 

x = —x + y 
ij = a(2A —1 )-y. 

Here, tangencies occur when y = x = 0, so there is no sliding region. Instead, as in Section 2.1, we 
have saddle behavior through two visible tangencies occurring at the same point, y = 0. 
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X 


(A) 



X 


(C) 


Figure 6. Example phase diagrams of (fTUD . with a = 1 and (a) k < 0, (b) k = 0, 
and (c) k > 0. As before, red curves represent the trajectories that have behavior 
similar to an unstable manifold, and blue curves exhibit behavior similar to a stable 
manifold. The points show the equilibria of the system, and the squares show where 
the sliding region begins/ends. Dashed black lines represent the nullclines of the 
systems. 


4.2. Energy balance model. In this section we analyze the energy balance model (JT]) from the 
introduction. We will reproduce the equations for reference: 

T =Q(l-a(T))-(A + BT) 

A = —mA + nT + C, 


where 


a{T) = 


+ or 


Oii QLq, 


tanh 


T- T 0 
D 


As we will see, the system has a repelling sliding region as D —» 0. The system is not quite in the 
same form as the system in Theorem 2 as is, however there is a change of variables by which it is 
possible. The resulting system is no more enlightening than © is in its current form, so we forgo 
the change of variables and analyze the system as is. Computing the Jacobian gives 


( 11 ) 


J = 


—Qa'(T) - B -1 
n —m 


and det(J) = m(Qa'(T ) + B) +n. The dependence of a on T is through a hyperbolic tangent with 
a negative coefficient (a, — a w )/2. In the limit as D —> 0, a(T) limits to 


a 0 (T) = | 


if T < T 0 
if T > T 0 . 


( 12 ) 
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Define Ai = Q( 1 — a*) — BTq and define A w similarly. If the A-nullcline intersects the line T = Tq 
between A{ and A w , then Figure [T] indicates there will be an equilibrium near To at say, T*(D), 
that limits to To as D — > 0. This occurs if 

Tl 

Ai < —To + C < A w . 
m 

Also, in the D = 0 limit, the derivative of a(T*(D)) —> —oo. So for D small enough, the equilibrium 
at T*(D) will be a saddle. As shown in Figure [71 ([Tj) limits to a system with a repelling sliding 
region with an attracting critical point. The two proper equilibria that remain in the nonsmooth 
system are fully attracting since at these points det(J) = mB + n > 0 and Tr( J) = —{B + m) < 0. 
Thus, the zombie saddle at Tq serves as an unstable object between the two attracting equilibria. 



Figure 7. Phase diagram of ([3]) in the nonsmooth limit. Parameter values are the 
same as in Figure [U i.e. m = n = 1, Q = 300, i = 0.8, w = 0.2, To = 290, B = 4, 
and C = 1300. Blue curves represent stable manifold-like behavior, red curves 
correspond to unstable manifold-like behavior, and the dashed lines represent the 
nullclines of the system. The points show the equilibria of the system, and the 
squares show where the sliding region begins/ends. 


5. Discussion 

It has been established that there are perils associated with using a nonsmooth approximation 
of a smooth dynamical system [Tl], and one might view this paper as yet another chapter of that 
cautionary tale—after all, a nonsnrooth system may not even reflect the number of equilibrium 
points of the system it is supposed to approximate. On the other hand, we view this as highlighting 
similarities between certain smooth and nonsmooth systems. In this case, there is saddle behavior 
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in the nonsmooth system that mimics what is observed in the smooth system. The saddle behavior 
is caused either by a double-tangency along the splitting line as in (J7|) or by a an equilibrium of a 
sliding solution as in systems (JHJ) and ([9]). We refer to the limiting points in the nonsmooth systems 
as having “zombie behavior” to reflect that even though the saddle equilibrium was destroyed, the 
saddle behavior still exists in an altered state. 

The recognition of the saddle behavior, especially the existence of trajectories analogous to 
stable manifolds of a zombie saddle, can be a useful tool for analyzing a dynamical system. Under 
certain parameters, the motivating conceptual climate model (JT]) exhibits bistability. In the smooth 
system, there is a saddle point between the two attracting equilibria, and the stable manifold of the 
saddle separates the basins of attraction of the stable equilibria. Analysis of this sort is standard 
in smooth systems. The existence of a zombie saddle allows us to carry out similar analysis in the 
corresponding nonsmooth system. Future exploration may examine the manner in which (un)stable 
manifold of a saddle in the smooth system limits to the (un)stable manifold of the zombie saddle. 
The computational difficulty in this endeavor arises because the zombie saddle is not truly an 
equilibrium, and thus dynamics near the point are not described by its linearization. Therefore, 
we cannot rely on the eigenvalues and eigenvectors of the linear system to provide the desired 
information. 

In smooth systems, if the stable manifold of a saddle is co-dimension one, then it forms a 
separatrix, dividing phase space into regions with different dynamics. For this reason, it would be 
useful to find similar results for zombie saddles in higher dimensions where computing basins of 
attraction is generally more difficult. 
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